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Abstract 

We  analyze  in  detail  conversion  of  a  pump  wave  into  two  daughter  waves  (DW’s)  on  an 
infinitely  thin  nonlinear  layer,  assuming  that  the  conversion  is  mi*’ated  by  weak  seed  pulses  of 
the  DW's.  Three  cases  corresponding  to  different  signs  of  the  group  velocities  of  the  DW’s  aure 
considered  separately.  In  the  case  when  both  velocities  are  opposite  to  that  of  the  pump  wave, 
a  threshold  value  of  the  ’’strength”  of  the  nonlinear  layer  is  found,  above  which  oscillations 
between  the  conversion  and  the  inverse  process  of  recombination  set  in.  In  the  case  when  the 
DW  group  velocities  have  different  signs,  the  shapes  of  the  generated  DW  pulses  do  not  depend 
upon  shapes  of  the  corresponding  "seed”  ptilses.  In  this  case,  multiple  solutions  are  possible, 
provided  the  nonlinearity  is  strong  enough. 

1  Introduction 

It  is  well  known  that  three-wave  resonant  interactions  (3WR1)  play  a  fundamentad  role  in  a  number 
of  problems  of  plasma  physics  and  nonlinear  optics  [1].  In  many  cases,  it  is  necessary  to  analyze 
these  interactions  in  inhomogeneous  media  [2,3].  An  important  physiczd  realization  of  the  3WRI  in 
an  inhomogeneous  medium  is  interaction  of  an  intensive  electromagnetic  wave  with  a  plasma.  The 
linear  interaction  of  electromagnetic  waves  with  an  inhomogeneous  medium  is  a  traditional  problem 
of  nonlinear  optics  [4]  and  plasma  physics  (see,  e.g.,  refs.  [5, 6,7, 8]).  Recently,  attention  has  been 
attracted  to  resonant  nonlinear  interactions  of  the  waves  with  a  strongly  inhomoneous  plasma  layer 
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(sheath),  which  can  find  a  promising  practical  application  to  design  means  of  communication  with 
a  space  vehicle  reentering  the  atmosphere  [9]. 

In  ref.  [3],  the  smalysis  WcU5  developed  for  two  situations:  When  the  inhomogeneity  was  very 
smooth,  and  when,  contrary  to  this,  the  nonlinearity  was  residing  in  a  very  narrow  layer,  so  that 
the  nonlinear  terms  in  the  governing  equations  were  proportional  to  a  delta-function.  The  objective 
of  the  present  paper  is  to  complete  the  analysis  for  the  thin  layer  started  in  ref.  [3]. 

The  system  of  equations  governing  the  3WRI  in  an  infinitely  thin  layer  has  the  form  [3] 


(t^i)<  +  vi(ui)*  =  -u<5(x)upu5. 

(1.1a) 

Mt  +  t^(t‘2)*  =  -i€S{x)UpU\, 

(1.16) 

(Up)t  +  v{ttp):i  =  -tc5(x)UiU2, 

(l.lc) 

where  ui{x,t),U2{x,t),  and  Up{x,i)  are  complex  envelopes  of  the,  respectively,  DW’s  and  pump, 
t;i,U2,  and  v  being  their  group  velocities.  In  what  follows,  the  pump  velocity  v  will  be  assumed 
always  positive  by  definition,  while  the  DW  velocities  wi,2  may  have  either  sign. 

In  an  infinitesimal  vicinity  of  the  point  x  =  0,  one  may  neglect  the  temporal  derivatives  in  eqs. 
(1),  reducing  them  to  the  system  of  ODE’s  with  singular  coeflScients.  Even  in  the  case  of  a  finitely 
thin  layer,  this  approximation  is  still  valid,  provided  the  envelope  is  slowly  varying  in  time,  such  that 
Ax|u(l  <<  Ivu]  where  Ax  is  the  width  of  the  thin  layer.  We  also  remark  that  if  we  interchange  x  and 
t,  we  see  that  this  problem  is  also  solvable  by  the  1ST  for  the  {3WRI)  [1].  Upon  this  interchange, 
we  now  have  the  initial  value  problem  where  the  profiles  are  specified  at  t  =  0",  they  then  interact 
for  a  short  time  with  a  very  large  coupling  constant  imtil  t  =  O'*',  at  which  time  we  wish  to  know 
the  final  profiles.  Although  the  1ST  can  formally  solve  this  problem  the  complexities  of  the  inverse 
scattering  equations  can  be  bypassed  with  the  following  simpler  analysis. 

An  important  part  of  this  problem  concerns  the  relative  phase  of  each  envelope^.  In  general, 
each  envelope  has  an  amplitude  and  a  phase, 


(1.2a) 
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Up  =  (l-2c) 

Then  from  (1.1),  upon  neglecting  the  time  derivatives,  we  have 


—  £5(i)i4pj42  sin  a 

(1.3a) 

V2-^  =  €5(x)>4pi4i  sin  a 
ax 

(1.36) 

=  — €(5(x)AiA2sina 
dx 

(1.3c) 

ui— =  —edix)—^ — COSQ 
dx  Ai 

(1.4a) 

doi2  r/  X  ApAi 

1)2—7—  =  —edlX)—^ — COSQ 

dx  A2 

(1.46) 

dcxp  .  A1A2 

W— 7-^  =  —€d{x) — ; COSQ 

dx 

(1.4c) 

where 


a  =  Op  —  Qi  —  Q2 


(1.5) 


Now  from  (1.4) 


4  ) 

A1A2 

V2A2/ 

vAp 

COSQ 


(1.6) 


Note  that  the  smallest  amplitude  in  (1.6)  will  determine  to  which  fixed  point  (a  =  ±|)  the  system 
will  be  driven  toward.  If  one  amplitude  is  much  smaller  than  the  other  two,  then  the  sign  of  its 
group  velocity  determines  which  fixed  point  is  stable  and  a  will  be  driven  very  rapidly  to  that  stable 
point.  If  the  two  smallest  amplitudes  are  comparible,  then  if  the  signs  of  their  group  velocities  are 
different,  they  can  compete  for  determining  the  stable  point.  In  any  case,  any  value  of  a  away 
from  ±1  is  unstable  and  it  is  only  near  the  values  of  ±|  that  each  of  the  individual  phases  become 
stationary  [see  (1.4)].  We  shall  therefore  assiune  that  a  =  and  which  value  we  use  shall  be  the 
same  as  the  stable  fixed  point  of  (1.6).  Of  course,  these  two  stable  points  only  differ  by  a  phaise  of 
TT,  which  simply  means  a  sign  chztnge  of  the  amplitude.  Sometimes  we  may  find  it  convenient  to 
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simply  let  the  amplitude  carry  the  sign  instead  of  worrying  about  the  phase  change,  particiilarly 
when  the  amplitude  pcisses  through  zero. 

As  it  immediately  follows  from  eqs.  (1-3),  one  can  eliminate  the  amplitudes  Ai  and  A2  in  favor 
of  A,  [4] 


(1-7) 

Cn{t)  being  arbitrairy  functions  of  time.  Actually,  eq.  (1.7)  is  nothing  but  the  well-known  Manly- 
Rowe  relation  [1].  At  last,  noting  that  Ai  and  A2  are  each  positive,  inserting  eq.  (1.7)  into  eq. 
(1.3c),  one  obtains  the  first-order  equation  for  Ap{x) 

^  =  -€S{x)suxay/(viV2)-^{Ci  -  Al)iC2  -  Aj)  (1.8) 

upon  noting  that  u  >  0. 

A  general  problem  that  we  are  going  to  discuss  is  to  use  eq.  (1.8)  to  analyze  the  conversion 
of  an  arbitrary  initial  pump  pulse  u^'^(t  —  f  )•  To  do  this,  it  is  necessary  to  specify  the  functions 
Cn(t)(n  =  1,2)  in  eq.  (1.8),  which  one  can  do  considering  boundziry  conditions  at  x  =  ±0.  In  turn, 
this  analysis  proves  to  be  essentially  difi’erent  for  different  signs  of  the  DW  velocities  Therefore, 
in  what  follows  we  will  analyze  separately  three  different  cases. 

2  Forward  Conversion:  Vi,V2>  0 

For  the  finite  slab  case,  this  case  of  forward  conversion  could  be  either  the  decay  CEise  (0  <  <  v  < 

V2  or  0  <  V2  <  <  ^i)  or  the  stimulated  backscatter  (SBS)  case  (0  <  u  <  i;i,t;2  or  0  <  t;i,t;2  <  u) 

depending  on  whether  the  magnitude  of  velocity  of  the  pump  lies  between  or  outside  of  the  velocities 
of  the  DWs.  However,  these  distinctions  and  differences  become  irrelevant  in  the  limit  of  a  thin 
layer,  with  the  decay  case  and  SBS  case  having  the  same  limit.  Here,  we  expect  that  the  initial 
pump  pulse  —  ^)  will  be  converted,  in  part,  into  DW  pulses  u„  = 

give  rise  to  a  residual  (final)  pulse  of  the  pump  wave  —  ^).  To  initiate  the  conversion  process, 
it  is  necessary  to  laxmch,  together  with  u^*^(t  —  ^),  small  "seed”  DW  pulses  —  f-)  (otherwise 
we  will  have  the  trivial  solution  =  0).  Now  for  the  value  of  q.  In  this  case,  since  we  initially 
only  have  small  seed  DW  pulses,  the  last  term  in  (1.6)  can  be  neglected.  Then  a  simple  analysis 
easily  shows  that  q  =  |  is  the  stable  value  and  (1.8)  gives  that  the  pump  will  initially  deplete.  As 
it  depletes,  the  amplitudes  can  pass  through  zero.  If  that  happens,  we  simply  set  a  =  |  and  allow 
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the  amplitudes,  if  they  pass  through  zero,  to  carry  a  sign.  Next,  the  functions  Cn{t)  can  be  found 
from  eqs.  (1.7)  taken  at  x  =  0“  (i.e.  before  the  conversion): 


C„(()  =  s  +  jUM,n  =  1.2,  (2.1) 

where  the  initial  amplitudes  ^4^'^  and  are  related  to  the  initial  pulses  of  the  pmnp  and  daughter 
waves  by  eqs.  (1.2). 

According  to  what  was  said  above,  we  assume  that 


Un  «  (AW)^  (2.2) 

Next  one  should  insert  eq.  (2.1)  into  eq.  (1.8)  and  integrate  the  latter  equation,  making  use 
of  the  relation  S{x)dx  =  1.  This  must  yield  the  profile  of  the  residucil  (final)  pulse 
determined  by  A^'^(t)  and  A[*2(0'  ^  ^  general  case,  this  expression  involves  elliptic  integrals  and  is 
rather  complicated  [6].  However,  the  integral  may  be  simplified  if  one  makes  use  of  the  inequality 
(2.2).  Eventually,  we  find 

(2,3) 

where 


,  ^  VlUlit)+V2U2{t) 

^  ^  “  8u(A<*)(t))2  ®*P 


2eAW(t) 


fViV2 


(2.4) 


Although  the  pre-exponential  factor  in  the  expression  (2.4)  is  small  according  to  eq.  (2.2),  its 
smallness  can  be  readily  compensated  by  the  exponentied  if  c  is  sufficiently  large,  and  the  conversion 
rate,  proportional  to  K{t),  can  thus  be  nonsmall.  Note  that  A^^^(t)  vanishes  when  K{t)  attains  the 
value  unity,  and  with  subsequent  increase  of  the  amplitude  A^^^  is  again  nonzero  with  the 
opposite  sign,  attaining  the  value  A^^\t)  —  —  A*’^(t)  at  A^*^(t)  =  oo.  Although  the  amplitudes  A 
have  been  defined  in  eqs.  (1.2)  as  positive  quantities,  its  negative  values  imply  that,  beyond  the  full 
conversion  point  (A*^^  =  0)  the  DW’s  produced  by  the  decay  of  the  pump  start  to  recombine  into 
it.  As  it  was  explained  in  ref.  [3],  one  can  actually  expect  oscillations  between  the  pump  and  DW’s, 
if  the  delta-function  layer  is  replaced  by  a  structured  profile  with  a  finite  width.  In  the  framework 
of  the  present  analysis,  the  oscillations  will  appear  if  one  lauches  finite  but  small  seed  DW  pulses 
instead  of  the  vanishingly  small  ones  considered  above.  In  this  case,  the  approximation  (2.4) 
is  no  longer  valid  and  would  have  to  be  replaced  by  the  appropriate  elliptic  integral  which  woiild 
contain  this  oscillatory  behavior. 
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With  the  quantity  given  by  eqs.  (2.3)  and  (2.4),  one  can  find  the  DW  pulses  gen¬ 

erated  by  the  piunp  conversion.  iFYom  eqs.  (2.1)  and  (1.7),  taken  at  x  =  O'*’  (after  the  conversion), 
it  follows  that 


(XWI(t))“  «  f  ((A«(t))»  -  (2.5) 

The  DW  pulses  will  propagate  in  the  region  x  >  0  in  the  form  u„  =  Un(t  —  ^).  Thus,  in  the 
case  when  Vi  and  V2  are  positive,  the  pump  pulse  passing  through  the  nonlinear  layer  gives  rise  to 
a  residual  pump  pulse  and  to  DW  pulses,  all  propagating  in  the  szune  direction. 

At  last,  it  is  relevant  to  ancdyze  if  the  phase  stability  condition  remains  valid  during  the  ampli¬ 
tude  conversion.  Inserting  eqs.  (1.7)  and  (2.1)  into  (1.6),  in  which  Un  is  neglected  in  comparison 
with  (A^'^)^,  one  can  readily  find  that  the  fixed  point  a  =  |  is  no  longer  stable  if 

(2.6) 

i.e.,  when  the  conversion  is  too  deep  (this,  in  p<uticular,  means  that  the  point  of  the  full  conversion, 
at  which  vanishes,  lies  in  the  unstable  domain).  Using  eqs.  (2.3)  and  (2.4),  it  is  easy  to  find 
in  an  explicit  form  a  boimdary  between  stable  and  unstable  domains  in  the  parametric  space.  In 
this  work  we  will  not  consider  how  the  solutions  should  be  modified  to  take  account  of  the  phase 
instability  in  the  domain  (2.6). 

3  Backw^ard  Conversion:  vi^V2  <0 

The  backward  conversion  case  corresponds  only  to  SBS  in  the  homogeneous  case  of  a  finite  slab. 
There,  the  decay  of  the  piunp  occurs  through  a  backscattering  process  whereby  the  pump  is  partially 
reflected  [1],  which  in  a  sense,  is  a  mode  conversion.  In  this  case,  the  generated  DW  pulses  will 
propagate  backwards,  i.e.  they  will  look  like  reflected  waves.  To  generate  them,  it  is  necessziry 
to  launch  small  "seed”  (initial)  pulses  at  x  =  O'*’,  Actually,  the  conversion  will  proceed  from 
X  =  0'*'  to  X  =  0".  With  regard  to  this,  eq.  (1.6)  with  negative  Vi  and  V2  tell  us  that  the  stable 
value  of  the  phase  a  is  always  Thus,  in  the  present  case  the  amplitude  evolution  equation  (1.8) 
also  gives  rise  to  depletion  of  the  inital  pump  pulse. 

According  to  what  was  said  above,  the  functions  C„(t)  are  determined  by  eq.  (1.7)  tzdcen  at 

X  =  0'*’; 
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(3.1) 


C„(t)  =  +  ^(AS'(f))"  s  +  jVM 

cf.  eq.  (2.1).  This  time,  it  is  assiuned 


t/„  «  (A<«)*  ,  (3.2) 

cf.  eq.  (2.2).  Note  that  with  this  approximation,  we  are  assuming  that  complete  conversion  will  not 
occur.  The  subsequent  analysis  can  be  developed  essentially  as  it  was  done  in  the  preceding  section. 
We  insert  eq.  (3.1)  into  eq.  (1.8),  and  then,  to  calculate  the  resultant  integral  approximately,  we 
make  use  of  the  inequality  (3.2).  Eventually,  we  arrive  at  the  following  trzmscendental  equation 
relating  to 


(|vi|t/i  +  h|C/,)(.4«  +  .4l/»)'  ’  ~ 


(3.3) 


Elementary  analysis  of  this  equation  demonstrates  that  there  is  a  threshold  value  Cthr  such  that 
eq.  (3.3)  has  two  solutions  at  e  <  €</,,,  and  no  solution  at  e  >  €thr-  One  can  find  an  approximate 
expression  for  Cthr  ouce  again  making  use  of  the  inequality  (3.2), 


j2  _  8vuit;2 

"  KI£/i  +  |t^|£/2- 

Note  that  the  threshold  value  given  by  eq.  (3.4)  does  not  depend  on  and  it  is  large  in 
virtue  of  the  smallness  of  Ui  and  U2-  At  c  >  €thr,  the  above-mentioned  oscillations  between  the 
pump’s  decay  into  the  DW’s  and  the  reverse  process  of  recombination  of  the  DW’s  into  the  pump 
take  place,  so  that  €thr  is,  as  a  matter  of  fact,  the  threshold  for  the  onset  of  oscillations.  However, 
the  oscillations  czinnot  be  considered  in  the  fi-amework  of  the  approximation  used.  To  do  this,  it  is 
necessary  to  consider  the  excict  elliptic  integral  corresponding  to  eq.  (1.8). 

The  existance  of  two  solutions  below  the  threshold  implies  that  one  solution  should  be  stable 
and  one  unstable.  The  stable  branch  corresponds  to  >  0,  and  the  unstable  one  to  <  0.  In 
particular,  in  the  region  >  3^^  the  stable  branch  can  be  found  in  the  following  approximate 
form: 


/ViV2 


In 


8i;(A<‘))2 


2c  \V\\Ui-\-\v2\U2 
In  this  region,  this  branch  lies  in  a  region  of  small  A^^K 


(3.5) 
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To  conclude  this  section,  let  us  emphasize  once  again  that,  in  the  case  considered,  the  value 
remauns  the  stable  fixed  point  of  eq.  (1.6)  at  all  the  values  of  the  parameters. 

4  Mixed  Conversion:  vi  >  0,  t;2  <  0 

This  case  is  another  example  where  the  interaction  could  be  the  soltion  decay  case  (vi  <  v)  or 
the  stimulated  backscatter  case  {y  <  Vi).  In  this  case  where  vj  and  V2  have  different  signs,  the 
seed  pulses  of  the  two  daughter  waves  should  be  taken,  respectively,  at  i  =  O'*'  and  at  x  =  0". 
Accordingly,  the  functions  C\  and  C2  are  determined  by  eqs.  (1.7)  taken  at  x  =  O'*’  and  at  x  =  0”. 
To  determine  stable  values  of  the  phase  a,  one  can  consider  the  partictilar  cases  when  either  of 
the  seed  components  is  nonzero,  while  the  other  one  is  absent.  If,  for  instance,  only  the  first 
component  (n  =  1)  is  present,  the  conversion  process  proceeds  from  x  =  0“  to  x  =  O'*',  and  eq. 
(1.6)  with  positive  Ui  defines  a  =  |  as  the  stable  fixed  value.  It  is  easy  to  see  that  the  Scime  is  true 
in  the  opposite  particular  case,  when  only  the  second  seed  component  is  different  from  zero.  Thus, 
as  well  as  in  the  two  preceding  cases,  eq.  (1.8)  gives  rise  to  depletion  of  the  initial  pump.  However, 
the  important  difference  from  the  previous  cases  is  that  this  time  the  elliptic  integral  following  from 
eq.  (1.8)  converges  even  when  the  initial  seed  values  Aj,’'  are  exactly  equcd  to  zero.  Actually,  this 
means  that  the  results  which  will  be  put  forward  below  are  valid  provided  that 

(^(i))2_(^(/))2»(^W)2^  (4.1) 

Neglecting  the  ’’seeds”  according  to  what  was  said  above,  we  obtain 

C,(l)  =  C,(t)  =  (4.2) 

cf.  eqs.  (2.1)  2uid  (3.1).  Then,  it  is  straightforward  to  insert  eqs.  (4.2)  into  eq.  (1.7)  to  obtain  from 
it  the  relation  between  and  A*  in  the  following  general  form: 

6(A«)2  =  (1  +  2N)y/\^\K{y/l  -  [A(/)/A(*)]2),  (4.3) 

K{q)  being  the  complete  elliptic  integral  of  the  first  order  with  the  modulus  q  =  yjl  — 

The  integer  N  =  0,1,2,...  takes  into  accoimt  explicitly  the  above-mentioned  circumstance  that 
the  pump  and  daughter  waves  may  accomplish  N  cycles  of  the  full  conversion  in  the  course  of  the 
interaction  in  the  nonlinear  layer.  It  will  be  demonstrated  below  that  N  cannot  take  indefinitely 
large  values;  nevertheless,  we  conclude  that,  in  the  case  when  N  may  take  more  than  one  vzilue,  the 
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three-wave  conversion  problem  has  multiple  solutions.  At  the  same  time,  one  should  remember  that 
each  inversion  of  the  pump  decay  into  the  opposite  process  of  the  piunp  recombination  introduces 
a  jump  Aa  =  tt  into  the  phase  a  defined  in  eq.  (1.5).  So,  stability  of  different  branches  of  the 
multiple  solution  remcuns  to  be  clarified. 

To  represent  the  final  amplitude  as  a  function  of  the  initial  one,  we  introduce  the  function 
q  =  q{K)  as  an  inverse  to  K  =  K{q).  Then,  the  representation  following  from  eq.  (4.3)  is 


=  (A<‘))^ 


V(H-2JV)Y^|t/iV2|/. 


(4.4) 


Note  that  eqs.  (4.3)  and  (4.4)  demonstrate  correct  asymptotic  behavior  in  the  limit  e  -*  oo  {at  N 
fixed):  (f{oo)  =  1,  so  that  =  0,  i.e.  the  full  conversion  takes  place  in  this  limit.  On  the  other 
hand,  eq.  (4.3)  has  no  solution  in  the  region 


<  5(1  +  2N).  (4.5) 

Thus,  for  ziny  given  e  we  obtain  a  finite  number  of  solutions  of  the  mixed  conversion  problem.  In 
particular,  there  is  no  solution  of  the  form  sought  for  if  inequality  (4.5)  holds  for  N  =  0.  Actually, 
the  basic  assumption  made  above,  i.e.  omission  of  the  "seeds”  is  not  legitimate  in  the  region 
(4.5),  so  that  a  solution  in  this  region  must  take  a  more  complicated  form.  One  can  also  understand 
this  from  the  viewpoint  of  the  soliton  decay  case  discussed  in  Ref.  [1].  The  solution  for  soliton 
decay  is  very  sensitive  to  the  "seeds”,  with  the  scalelength  for  the  decay  linearly  proportionally  to 
the  natural  logeirithm  of  the  seed  amplitude.  Also  the  relative  phases  of  the  decay  products  can 
be  strongly  affected  by  the  relative  phases  and  amplitudes  of  the  seeds.  Consequently  it  is  not 
surprising  that  if  we  take  the  limit  of  vanishing  seeds,  we  fiind  a  multitude  of  possible  solutions.  We 
also  note  that  in  the  SBS  case  [1],  that  problem  is  a  two-point  boimdary  value  problem  as  is  here. 
Such  problems  can  be  remarkable  sensitive  to  initial  conditions,  particularly  near  any  homogeneous 
solutions.  That  may  also  provide  a  means  of  understanding  these  multiple  soltuions. 

At  last,  the  final  form  of  the  DW’s  produced  by  the  conversion  of  the  piunp  can  be  obtained 
from  eqs.  (4.2).  The  waves  ui  and  U2  propagate,  respectively,  as  transmitted  and  reflected 
waves  generated  by  the  nonlinear  layer.  The  important  fact  is  that  outside  of  the  region  (4.5),  the 
shapes  of  the  generated  DW  pulses  do  not  depend  on  the  form  of  the  "seed”  pulses. 

Coming  back  to  the  phase  stability  problem,  one  can  readily  find  that,  apart  from  the  phase 
jumps  at  the  decay-recombination  inversion  points,  the  stability  exponent  determined  by  eq.  (1.6) 


e(A<»>)^ 

\/|viV2l 
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changes  its  sign  at  the  point 


I 


+  ^(AW)4  +  (^(/))4  _  (>lC))2(^(/))2]  .  (4.6) 

u 

It  is  easy  to  check  that  the  criticai  v-\iae  of  A^  determined  by  eq.  (4.6)  always  belongs  to  the  interval 
(4*^^)^  <  .4^  <  (i4^’^)^,  so  that  «his  value  is  automatically  attained  in  the  course  of  the  conversion. 
Thus,  in  any  case  the  conversion  index  given  by  eq.  (4.4)  pertains,  strictly  speaking,  to  a  solution 
which  is  not  globally  stable.  Nevertheless,  this  result  may  be  quite  useful,  as  it  should  be  relevant 
at  least  in  tb**  case  when  the  initial  rate  of  the  phase  disturbances  is  low  enough. 
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